Algorithms for rank-de cient and ill-conditioned Toeplitz least-squares and QR factorization
نویسندگان
چکیده
In this paper we present two algorithms-one to compute the QR factorization of nearly rank-deecient Toeplitz and block Toeplitz matrices and the other to compute the solution of a severely ill-conditioned Toeplitz least-squares problem. The rst algorithm is based on adapting the generalized Schur algorithm to Cauchy-like matrices and has some rank-revealing capability. The other algorithm is based on adapting the augmented systems method to Toeplitz matrices. This algorithm does not suuer from the numerical inaccuracy of most Levinson-and Schur-based schemes proposed in the literature thus far because it avoids forming the normal equations either implicitly or explicitly. Some comparisons with more recent Toeplitz least-squares algorithms are also presented.
منابع مشابه
Parallel Algorithms for Toeplitz Systems
We describe some parallel algorithms for the solution of Toeplitz linear systems and Toeplitz least squares problems. First we consider the parallel implementation of the Bareiss algorithm (which is based on the classical Schur algorithm). The alternative Levinson algorithm is less suited to parallel implementation because it involves inner products. The Bareiss algorithm computes the LU factor...
متن کاملA weakly stable algorithm for general Toeplitz systems
We show that a fast algorithm for the QR factorization of a Toeplitz or Hankel matrix A is weakly stable in the sense that RR is close to AA. Thus, when the algorithm is used to solve the semi-normal equations RRx = A b, we obtain a weakly stable method for the solution of a nonsingular Toeplitz or Hankel linear systemAx = b. The algorithm also applies to the solution of the full-rank Toeplitz ...
متن کاملSuperfast and Stable Structured Solvers for Toeplitz Least Squares via Randomized Sampling
We present some superfast (O((m + n) log(m + n)) complexity) and stable structured direct solvers for m× n Toeplitz least squares problems. Based on the displacement equation, a Toeplitz matrix T is first transformed into a Cauchy-like matrix C, which can be shown to have small off-diagonal numerical ranks when the diagonal blocks are rectangular. We generalize standard hierarchically semisepar...
متن کاملNew Fast Algorithms for Structured Linear Least Squares Problems
We present new fast algorithms for solving the Toeplitz and the Toeplitz-plus-Hankel least squares problems. These algorithms are based on a new fast algorithm for solving the Cauchy-like least squares problem. We perform an error analysis and provide conditions under which these algorithms are numerically stable. We also develop implementation techniques that signiicantly reduce the execution ...
متن کاملFast Qr Factorization of Low-rank Changes of Vandermonde-like Matrices
VANDERMONDE-LIKE MATRICES LUCA GEMIGNANI Abstract. This paper is concerned with the solution of linear systems with coe cient matrices which are Vandermonde-like matrices modi ed by adding low-rank corrections. Hereafter we refer to these matrices as to modi ed Vandermonde-like matrices. The solution of modi ed Vandermondelike linear systems arises in the approximation theory both when we use R...
متن کامل